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(j)  =  Optical  phase 

A(|)  =  Change  in  optical  phase 

^ip^i2  “  Pockel  coefficients 


Cylinder  and  Fluid 

Angle  relative  to  radial  coordinate  r 
Angular  frequency 


=  Axial  wavenumber  (27r/L) 

=  Axial  wavelength 

=  j={ 

=  Circumferential  order  number 

=  Radial  displacement  at  the  cylinder/fluid 
interface 

=  Bessel  function  of  the  first  kind 
=  Bessel  function  of  the  second  kind 
=  Modified  Bessel  function 

=  Modified  Bessel  function 

=  Hankel  function  of  the  first  kind 

=  Calculation  radius 
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A  THREE-DIMENSIONAL  DYNAMIC  ELASTICITY  SOLUTION  FOR  WAVE 
PROPAGATION  IN  A  TWO-LAYERED  INFINITE  VISCOELASTIC  SOLID 
CYLINDER  WITH  OUTER  FLUID  LOADING 

INTRODUCTION 

This  research,  which  is  a  precursor  to  the  case  of  a  two-layered  cylinder  with  inner  and  outer 
fluid  loading,  treats  the  two-layered  solid  cylinder  subject  to  forced  harmonic  vibration  within  the 
scope  of  the  linearized  theory  of  dynamic  elasticity.  In  previous  efforts,  the  statically  loaded 
solid  cylinder  was  documented  by  numerous  authors,  and  the  forced  vibration  response  of  the 
solid  cylinder  was  investigated  by  Holden,^  Kaul  and  McCoy,^  and  others. 

In  the  theoretical  development  presented  in  this  report,  the  case  of  nonaxisymmetric  vibration 
is  discussed  in  terms  of  circumferential  order  number  n.  Results  are  presented  for  the  first  two- 
order  numbers,  specifically  n  =  0  and  n  =  1.  Three  excitation  states  are  considered:  the  normal 
pressure  P^,  the  longitudinal  shear  stress  P^,  and  the  circumferential  shear  stress  Pq. 

The  development  proceeds  using  scalar  and  vector  displacement  potentials  to  describe  the 
motion  of  the  solid.  The  boundary  conditions  are  specified  in  terms  of  the  stresses  and 
displacements,  and  thus  must  be  written  in  terms  of  the  displacement  potentials.  For  the  case  of 
the  two-layered  cylinder,  the  problem  is  reduced  to  solving  a  lO-by-10  system  matrix  for  the 
undetermined  complex  coefficients.  This  matrix  requires  inversion  at  each  location  in  the 
wavenumber  and  frequency  plane,  allowing  all  dynamic  quantities  (stress,  strain,  and 
displacement)  to  be  calculated  from  the  coefficients.  Even  though  the  solution  is  closed  form,  a 
sizeable  numerical  computation  is  still  necessary  to  extract  information  from  the  equations 
governing  the  solution. 

An  interesting  application  for  the  models  derived  here  is  in  the  field  of  fiber  optics.  The 
response  of  an  optical  fiber  to  external  pressure  loading  has  been  the  subject  of  much 
investigation  in  the  literature.  For  example,  in  the  early  eighties  Lagakos^  developed  static 
elasticity  descriptions  for  the  behavior  of  a  coated  optical  fiber  subject  to  external  pressure 
loading.  More  recently,  Hull'*  developed  a  model  for  the  optical  phase  sensitivity  of  a  bare 
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optical  fiber  subject  to  an  axisymmetric  pressure  excitation  at  a  definite  wavenumber  and 
frequency.  What  has  been  lacking  in  the  field  of  fiber  optics,  however,  is  a  dynamic  model  for 
the  viscoelastic  coated  fiber  subject  to  both  axisymmetric  and  nonaxisymmetric  forced  vibration 
at  definite  wavenumber  and  frequency.  This  report  describes  such  a  model  and  presents  results 
for  a  cylinder  with  material  properties  and  dimensions  consistent  with  those  of  an  optical  fiber. 
Numerical  results  for  both  the  coated  fiber  (with  an  outer  diameter  of  8  mm)  and  bare  fiber  are 
presented. 
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THE  PHYSICAL  MODEL 

The  physical  model  under  analysis  is  shown  in  figure  1.  A  solid  cylinder  (cylinder  1)  of 
radius  a  is  joined  to  a  second  cylinder  (cylinder  2).  Cylinder  2  is  in  intimate  contact  with  cylinder 
1  at  radius  r  =  a,  extending  to  radius  r  =  b.  At  radial  distance  b,  cylinder  2  is  in  contact  with  the 
outer  fluid.  The  outer  fluid  extends  to  infinity  in  radial  coordinate  r.  The  cylinders  and  fluid 
extend  to  plus  and  minus  infinity  in  the  longitudinal  x-direction. 

As  stated  earlier,  the  system  is  treated  within  the  context  of  the  linearized  theory  of  elasticity. 
The  assumption  of  linearity  implies  that  both  the  cylinders  and  fluid  obey  linear  constitutive 
equations.  For  the  cylinders,  this  means  that  there  is  a  linear  relationship  between  stress  and 
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Strain.  For  the  fluid,  deviations  from  the  hydrostatic  pressure  are  a  linear  function  of  the 
condensation  (fractional  density  change)  and  similarly  the  dilatation  (volume  strain).  Strains  are 
assumed  to  be  small  compared  to  unity  for  both  the  cylinder  and  fluid;  the  response  therefore  is 
linearly  related  to  the  excitation.  Both  media  are  restricted  to  small  signal  analysis,  i.e., 
excitation  levels  on  the  acoustic  scale. 

Damping  is  included  by  the  use  of  a  complex  modulus  of  elasticity — ^the  well-known 
structural  damping  form.  Both  cylinders  are  treated  distinctly  because  each  possesses  different 
material  properties.  The  outer  fluid  is  modeled  as  inviscid  for  the  purpose  of  fluid  loading  the 
moving  structure.  The  radial  velocity  of  the  fluid  and  the  second  cylinder  are  identically  equal  at 
r  =  b.  Under  the  inviscid  assumption,  the  outer  fluid  interacts  with  only  the  radial  motion  of  the 
cylinder  and  is  uncoupled  from  the  circumferential  and  longitudinal  motions.  A  perfect  slip 
condition  is  allowed  between  the  cylinder  and  fluid  for  the  longitudinal  and  circumferential 
directions.  These  assumptions  avoid  boundary  layer  physics  within  the  context  of  this 
development. 

Three  harmonic  excitations  are  applied  by  the  outer  fluid  to  cylinder  2  at  r  =  ft:  a  radial 
pressure  P^,  shear  stress  P^,  and  shear  stress  Pq.  In  practice,  the  shear  stresses  P and  P q  would 
only  be  developed  by  the  viscosity  of  a  real  outer  fluid.  Therefore,  the  amplitude  distribution  (at 
r  =  b)  of  the  excitations  in  the  wavenumber  and  frequency  domain  must  be  known  a  priori', 
otherwise  a  calculated  quantity  is  normalized  by  an  excitation  stress  magnitude.  As  shown  in 
figure  1,  the  excitations  are  at  a  definite  wavenumber  {k)  and  frequency  (to)  and  are  of 
circumferential  order  number  n.  The  response  of  the  cylinder  will  therefore  be  a  function  of  k, 
to,  and  n,  as  well  as  of  the  geometric  and  material  properties  of  the  composite  system. 
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EQUATIONS  OF  MOTION 
DISPLACEMENT  POTENTIALS 

The  equations  governing  the  motion  of  a  homogeneous,  isotropic,  linearly  elastic  body  are 
formulated  with  the  strain-displacement  relations  (equation  (1)),  the  constitutive  relations 
(equation  (2)),  and  the  stress  equations  of  motion  (equation  (3))  of  reference  5: 


=  \%S,7  +  2(t^ey’  (2) 

and 

=  PA  ’  (3) 

where  Ey  are  the  components  of  the  strain  tensor,  Xy  are  the  components  of  the  stress  tensor,  m,  are 
the  components  of  displacement,  and  p,,  is  the  density  of  the  cylinder.  The  components  of 
acceleration  are  indicated  by  u-  and  the  body  forces  are  indicated  by  fj.  The  comma  denotes 
partial  differentiation  in  equations  (1)  and  (3)  as 


dUi 

dxj 


dxj 


(4) 


where  the  components  of  the  vector  of  spatial  coordinates  Xj  are  equal  to  Xj  =  r,X2  =  0,  and  x^  =  x. 

To  obtain  an  equation  only  in  terms  of  the  displacements  and  their  derivatives,  equations  (1) 
and  (2)  are  substituted  into  equation  (3),  resulting  in^ 


A  +  =  PA- 


(5) 


Equation  (5)  is  expressed  in  vector  notation  as 


|i  V  a  -H  (?l  +  p,  )  VV  •  «  =  p  u . 


(6) 
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The  displacement  vector  is  decomposed  using  the  scalar  and  vector  relationship 

M  =  Vtp  +  Vx\i/  ,  C7) 

where  cp  =  cp  (r,  0,  x,  t)  is  the  scalar  displacement  potential  and  \|/  =  \|/  (r,  0,  x,  r)  is  the  vector 
displacement  potential.  Equations  (6)  and  (7)  reduce  to 

V  ( +  2[i^)  V2(p  -  p^ip)  +  Vx  -  p^v|/)  =  0 .  (8) 

From  equation  (8),  we  obtain  the  following  scalar  and  vector  displacement  potentials. 

V2(p  =  \(p 
and 

=  4(f.. 

c, 

The  dilatational  and  shear  wave  phase  velocities  are  given  by 

2  +  2p.^  2  _  l^c 

Pc  Pc 

EcVc 

where  the  Lame  constants  and  are  equal  to  X^  =  ^ ^  ^ 

modulus  of  the  material,  respectively.’ 

For  the  present  cylindrical  geometry,  we  have  the  following  scalar  and  uncoupled  vector 


potentials: 

V2(p  =  ^(p 

(12) 

and 

=  -2V.- 

(13) 

(9) 

(10) 

(11) 

-  and  the  shear 

-2v,) 
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Two  coupled  vector  potentials  in  r  and  0  are 


^2  2 


2 

1^  M^e 

2^  2  ’ 
c^dt 


The  Laplacian  is  defined  in  cylindrical  coordinates  as 


y2  _  1  a  ,  1  9^ ,  8^ 

3r"  '■3'-  r"3e'  37 


Considering  the  scalar  potential  given  by  equation  (12),  we  will  assume  a  separation-of-variables- 
type  solution  in  r  and  0.  The  scalar  potential  tp  (r,  0,  a:,  t)  is  then  of  the  form 


(p  =  Q(r)0(0)e 


i  {kx-m) 


Equation  (12)  yields  the  following  two  equations: 


d  Q  IdQ 


2  \  2 
--k  Q--Q 


2 

d  ©  2_  ^ 

— j  +  n  0  =  0. 

dQ 


The  solution  of  equation  (18)  is  in  terms  of  Bessel  functions.  For  a  solid  cylinder,  such  as 
cylinder  1,  the  field  quantities  must  be  finite  at  the  center  (r  =  0);  therefore  we  retain  only  Bessel 
functions  of  the  first  kind,  J^.  However,  for  a  hollow  cylinder,  such  as  cylinder  2,  Bessel 
functions  of  both  the  first  and  second  kind  are  retained.  The  solution  to  equation  (19)  is  in  terms 
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of  the  sine  and  cosine  of  argument  nG.  The  total  solution  is  therefore 

,  T  T  /  X  ^  (kX-COt) 

(p  =  [Ajcos  (n0)  +  A2Sin  (nG)  ]/^  (pr)  e  , 

with  p  defined  as 


A  similar  solution  for  the  x-component  of  the  vector  potential,  is 

^  /  X  ^  (lcx-(Ot) 

=  [B^cos{nQ) +82^1x1  {nQ)]J^(qr)e 

with  q  defined  as 


(20) 


(21) 


(22) 


(23) 


The  potentials  for  ¥e  can  be  uncoupled  by  noting  that  V|/^  and  \pe  contain  trigonometric 

functions  of  0  and  that  a  sine  dependence  on  0  in  ti/;.  is  equivalent  with  a  cosine  dependence  on  0 
in  Ve,  and  vice  versa.*  With  the  use  of  the  condition  V  •  \|/  =  0  to  eliminate  an  extra  constant, 
\j/^  and  xpe  become 
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Displacement  Potential  Summary  for  Cylinder  1 

In  summary,  the  displacement  potentials  for  the  solid  cylinder  are  given  in  terms  of  Bessel 
functions  of  the  first  kind,^  with  argument  p  being  based  on  the  dilatational  wave  phase  velocity 
and  argument  q  on  the  shear  wave  phase  velocity  as  follows: 

(p  =  ipr)  cos  (nS) 

V.  =  sin(n0)e'^““^\ 

=  -Ci/„^i(9r)cos(«0)e'^"^"“'\ 


Displacement  Potential  Summary  for  Cylinder  2 

The  outer  cylinder  requires  Bessel  functions  of  the  first  and  second  kind.  The  summary  of  the 
scalar  and  vector  displacement  potentials  is 

9  =  [A^J^ipr)  +A2F„(pr)]cos(n0)c'^^^““'\ 

V;-  =  lC,J„^,{qr)  +C2T^^i(^/-)]sin(n0)e''^""~“^ 

Ve  =  [-Ci/„  +  i(^r)-C2T„^i(<?r)]cos(n0)e'^^^““'\ 
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PARTICLE  DISPLACEMENT/DISPLACEMENT  POTENTIAL 

Intermediate  relationships  between  the  particle  displacements  and  the  displacement  potentials 
are  given  by  equations  (28)  through  (30),^®  where  is  the  displacement  in  the  r-direction,  is 
the  displacement  in  the  circumferential  0-direction,  and  is  the  displacement  in  the  longitudinal 


x-direction: 


dr  ^  r30  dx 


^  1^ 

^ r90^3x  dr 


8(0  1  a ,  , 

^  dx  rdr  ®  rdQ 


STRESS/DISPLACEMENT  EQUATIONS 

In  the  problem  posed,  the  boundary  conditions  are  prescribed  with  the  stresses  and 
displacements.  Equations  (31)  through  (36)  relate  the  stress  components  to  the  particle 
displacements.  We  now  have  a  complete  progression  from  displacement  potentials  to  stresses, 
which  will  allow  us  to  solve  for  the  motion  of  the  composite  system: 


rr  cl  y.  ^30  0^  I 


.  “c  .  .  O.. 


”  c  ,  ci  C,9,,  C  i  C 

8«  -  M  57  r  rae  3;t  J “^4  r  rSG  )' 


- 1 

f.  ^30  J 
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dr  r  rdQ 


(34) 


l8w  8v 

7ae 


(35) 


and 


du^ 

dx  dr 


(36) 


The  indicial  notation  is  interpreted  in  the  following  manner.  The  first  index  refrt-s  to  the  plane 
normal  to  the  coordinate  identified  by  the  first  index.  The  second  index  refers  to  the  direction  in 
which  the  stress  acts.  Stresses  t„.,  Tqq,  and  are  normal  stresses,  whereas  x,.^  is  a  stress  acting 
in  the  plane  perdendicular  to  the  radial  coordinate  in  the  radial  direction.  The  remaining  three 
stresses  (x^-e,  Xe;,.,  and  are  shear  stresses.  Shear  stress  x,.e  indicates  a  shear  stress  acting  in  the 
plane  perpendicular  to  the  radial  coordinate  in  the  circumferential  9-direction. 


STRAIN/DISPLACEMENT  EQUATIONS 

These  equations  will  be  required  for  the  particular  output  needed  in  the  simulations  that 
follow: 

du 

e..  =  ^  ’  (37) 


and 


^ee  “ 


r  r30 


(38) 


(39) 
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STRESS,  STRAIN,  AND  DISPLACEMENTS  IN  TERMS  OF 
DISPLACEMENT  POTENTIAL  EQUATIONS 


The  stress  equations  needed  to  compose  a  suitable  set  for  the  application  of  the  boundary 
conditions,  together  with  the  displacement  equations,  are  written  in  terms  of  the  displacement 
potentials  by  substituting  equations  (28)  through  (30)  into  equations  (31),  (34),  and  (36)  to  obtain 
equations  (40),  (41)  and  (42): 


rr 


=  x. 


/  2  2  2 


'  C  ±  '  T  c 

+  -—  + 


+  • 


3^2  rdr  r'^dS^ 


+  2a 


^  idxi/^  3v|/e 


+ 


3^2  ^230  rdQdr  dxdr 


(40) 


2  2  2  2  ^ 

^  rdQdr  ^  dxdr  rdx  rdxdQ  dr^  rdr  r^dQ^ 


(41) 


/  2  2  2  2  2 

5  9c  i5¥c  l5¥r  5  \i/0  5¥e  1  l^Ve 

3v2  ^  0r2  f2^0  rdr  rdOdx 


'c  I  -  rr 
■  ^ 


drdx  r^dQ  rdQdr  3x2  3^2  ^2 


(42) 


This  same  substitution  is  performed  on  strain  equations  (37),  (38),  and  (39)  to  obtain  equations 
(43),  (44),  and  (45): 


2  22 

l5  5  ¥9 

^rr  ~  ^230  rdQdr  dxdr  ’ 

2  2  2 

_  i59c  i5  9c  ^  1^  ¥c  l5¥;, 

^ee  ”  737  ''’72^  ~  rdx  ^  r^dQ'^  ^  rdxdQ  rdQdr' 

2  2  2 

5  9c  .  5  ¥e  ,  l5¥e  1^  ¥, 

^  ^  dxdr  ^  rdx  rdxdQ  ' 


(43) 


(44) 


(45) 
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Equations  (46),  (47),  and  (48)  are  repeated  to  aid  the  reader: 

dr  ^  rdQ  dx 


(46) 


rdQ  ^  dx  dr 


(47) 


and 

dtp,  a\|/9  1  ia\|/, 

dx  dr  raO 


(48) 
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OUTER  FLUID 

As  the  surface  of  the  cylinder  is  set  into  motion  by  the  excitation  of  a  pressure  wave  of 
magnitude  P^,  a  pressure  field,  p^,  is  generated  in  the  outer  fluid  by  the  cylinder  surface.  When 
the  longitudinal  wavenumber  associated  with  the  vibration  of  the  cylinder  surface  is  smaller 
than  co/c^,  the  pressure  radiated  into  the  outer  fluid  is  in  the  form  of  a  propagating  wave.  If  the 
longitudinal  wavenumber  of  the  cylinder  vibration  is  greater  than  (o/c^,  the  pressure  field  decays 
exponentially  from  the  surface  of  the  cylinder  in  the  positive  radial  direction.  The  fluid  is 
modeled  as  an  ideal  (nonviscous)  linearly  elastic  media  that  cannot  sustain  shear  stress,  even 
when  it  is  in  motion.^^  The  equation  of  state  is  then 


0,1  = 


whose  components  are  given  by  Hooke’s  law,  o^j  =  +  2|Jye^^.,  by  setting  |j^=  0: 


^11  ~  ^22  “  ^33  -  '^f^kk  ~  Pf 


^12  ~  ^13  “  *^23  ”  ®  ■ 


The  following  scalar  displacement  potential  describes  the  dynamic  motion  of  the  fluid  in 
cylindrical  coordinates:^^ 

=  0.  (51) 


where 


= 


ar"^^^"''r2a02‘'az2 


In  the  case  of  fluids,  cy  =  which  is  equivalent  to  where  P  is  the  bulk  modulus  of  the 

fluid.  The  particle  displacement  and  the  pressure  in  the  fluid  can  be  expressed  in  terms  of  the 
displacement  potential  as  follows: 


Uf  =  Vcp. 
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2 

d  (p. 

Pf  ~ 

^  dt 


For  propagating  wavenumbers,  the  scalar  displacement  potential  is  given 


(p^  =  (§2^  cos  ■?2 


CO  ,2 


For  an  outgoing  wave  described  by  e  ’  ,  the  field  in  the  outer  fluid  is  physically  a  radiation 

from  the  surface  of  the  moving  cylinder.  The  radial  dependence  in  equation  (55)  is  a  Hankel 
function  of  the  first  kind,  indicated  by  superscript  (1);  i.e.,  {g^r)  =  J„  (^2^)  • 


Differentiating  equation  (55)  with  respect  to  radius  and  time  yields  equations  (56)  through 
(58)  for  the  fluid  pressure,  particle  displacement,  and  particle  acceleration,  respectively: 


p  =  to^p^M//^^  (§2^)  cos  (nfl) 


u  =  {g^r)  cos 


/  =  (gjr)  cos  («0)  e 

)  n  ^ 


For  nonpropagating  wavenumbers  in  the  outer  fluid,  the  evanescently  decaying  field  is  given  by 
the  following  displacement  potential,  where  the  radial  dependence  is  in  the  form  of  the  modified 
Bessel  function  K‘ 


cp^  =  HK^  (f^r)  cos  (nG)  ^ 


i  (kx-dit) 


f2 


2  CO  for 


,2  CO 

k  >- 
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Similar  to  the  above  discussion  of  the  propagating  case,  the  pressure,  particle  displacement,  and 
particle  acceleration  are  given  by  equations  (60)  through  (62)  as 

if^r)  cos  (nO)  e  ,  (60) 

and 

2 

(f^r)  cos  («0)  .  (62) 

Figure  2  connects  the  wavenumber-frequency  plane  with  the  time-space  characteristics  of  the 
field  properties  in  the  outer  fluid.  The  vertical  axis  represents  wavenumber  with  respect  to  the 
longitudinal  x-axis  of  the  cylinder.  The  horizontal  axis  is  the  frequency  axis.  The  positive 
wavenumber  direction  corresponds  to  propagation  in  the  positive  x-direction  and  the  negative 
wavenumber  direction  corresponds  to  propagation  in  the  negative  x-direction.  There  are  two 
distinct  regions  in  the  wavenumber-frequency  plane  corresponding  to  propagating  and  non¬ 
propagating  fields  in  the  outer  fluid. 

To  review,  variable  k  refers  to  the  wavenumber  that  exists  on  the  surface  of  the  cylinder.  The 
vibrating  cylinder  surface  will  initiate  a  propagating  field  in  the  outer  fluid  when  0^/c/  >  k^. 
When  this  condition  exists,  the  radial  displacement  (m^)  of  the  cylinder  surface  must  project  (m^’) 
onto  the  propagation  direction  of  the  field  in  the  outer  fluid.  For  a  given  frequency,  the 
wavelength  of  cylinder  vibration  must  be  greater  than  the  corresponding  wavelength  of  vibration 
that  the  fluid  will  propagate.  The  axisymmetric  n  =  0  vibration  produces  a  symmetric  field  in  the 
outer  fluid  as  a  function  of  0.  The  spatial  representation  of  a  point  ki  cOi,  in  the  region  of 
propagating  wavenumber-frequency  pairs  given  by  >  0 ,  is  shown  in  figure  3,  where  Lj  is 

the  wavelength  and  Cj  is  the  velocity  of  propagation.  In  the  limit  of  >  0 ,  the  radiated  field  is 
physically  a  series  of  cylindrical,  radially  spreading  wavefronts,  emanating  outward  from  and 
parallel  to  the  surface  of  the  cylinder.  The  intermediate  range  shown  in  figure  3  results  in  the 
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Figure  2.  Outer  Fluid  Propagating  and  Nonpropagating  Regions  of  the  Wavenumber- 

Frequency  Plane 
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wavefronts  propagating  out  along  conical  planes  forming  an  angle  a  with  the  longitudinal  v-axis 
of  the  cylinder  (figure  3).  The  direction,  a,  of  the  radiated  field  relative  to  the  longitudinal  A:-axis 
of  the  cylinder  will  be  given  by  equation  (63): 


a  =  acos—  =  acos—  =  acos  — 

Li  Cj 


Propagation  ceases  as  =  k^.  The  fluid  will  only  support  propagation  at  k^,  remembering  that 


Figure  3.  Relationship  Between  Cylinder  Vibration  and  Radiated  Pressure  Pg 


The  condition  k^  =  k^  requires  that  Cj  =  c^;  therefore  the  propagation  in  the  fluid  would  be 
along  the  longitudinal  x-direction.  However,  there  is  only  coupling  between  the  cylinder  motion 
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and  a  propagating  fluid  pressure  field  when  the  radial  displacement  u  can  be  projected  onto  the 
propagation  direction  of  the  fluid.  This  is  not  possible  when  ~  ^  ^ 

orthogonal  to  each  other;  this  condition  is  the  onset  of  the  nonpropagating  region  shown  in  figure 

2. 


The  Hankel  function  variation  of  the  outer  fluid  displacement  potential  can  be  observed  in  the 
inset  TIME-SPACE  diagram  in  figure  2.  This  diagram  depicts  the  radial  spatial  variation  of  the 
field  quantities  at  a  frozen  instant  of  time. 

As  a  line  of  constant  frequency  is  traversed  in  the  figure,  the  cylinder  wavelength  L]  becomes 
shorter,  ultimately  reaching  a  length  less  than  the  length  of  the  corresponding  fluid  wavelength 
that  supports  propagation.  Under  this  condition,  there  is  no  longer  a  component  of  the  radial 
displacement  of  the  cylinder  that  can  be  projected  onto  the  propagation  direction  of  the  field  in  the 
outer  fluid;  equation  (63)  is  no  longer  valid;  and  the  field  in  the  outer  fluid,  induced  by  cylinder 
vibration  ceases  to  propagate.  The  radial  variation  of  the  field  quantities  is,  according  to  the 
modified  Bessel  function  depicted  in  figure  2,  for  the  nonpropagating  regions.  All  field 
quantities,  such  as  the  displacement  potential  and  all  parameters  related  to  it,  undergo  the  same 
spatial  variation. 
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SOLID  CYLINDER/FLUID  INTERFACE 
BOUNDARY  CONDITIONS 

From  this  point  on  in  the  derivation  that  follows,  subscript  c  (Lame  constants)  will  be  replaced 
with  1  or  2  to  indicate  cylinder  1  or  2.  Superscript  Cl  or  C2  will  be  added  to  all  other  pertinent 

quantities  to  distinguish  the  first  from  the  second  cylinder.  The  solution  to  the  solid  cylinder  in 

Cl  Cl  Cl 

contact  with  an  outer  fluid  is  obtained  by  solving  for  the  constants  Aj  ,  and  Cj  in 
equations  (26)  and  (27)  and  for  M  in  equation  (55)  (or  H  in  equation  (59))  for  the  pertinent 
boundary  conditions  by  considering  the  state  of  stress  and  displacement  at  r  =  a  at  the  cylinder 
and  fluid  interface.  The  four  boundary  conditions  are  given  by  equations  (65)  through  (68)  as 
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The  four  boundary  conditions  at  r  =  (3  are  represented  diagrammatically  in  figure  4.  All 
boundary  conditions  are  in  effect  over  the  entire  area  of  contact  between  the  cylinder  and  the 
outer  fiuid. 


Figure  4.  Boundary  Conditions  for  the  Single  Solid  Infinite  Cylinder  and  Fluid 
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SYSTEM  MATRIX 

The  four  boundary  conditions  in  equations  (65)  through  (68)  are  assembled  in  matrix  form  in 
equation  (70); 


5/?(l,l)  SR  (1,2)  SR  (1,3)  SR  (1,4) 

SR  (2,1)  SR  (2,  2)  SR  (2,  3)  0 

SR  (3,1)  SR  (3,  2)  SR  (3,  3)  0 

SR  (4,  1)  SR  (4,  2)  SR  (4,  3)  SR  (4,  4) 

In  numerical  order,  each  boundary  condition  corresponds  to  the  row  number  of  system  matrix 
SR.  The  elements  of  the  SR  matrix  are  shown  in  equations  (71)  through  (86).  Solution  of 
equation  (70)  at  each  point  in  the  wavenumber-frequency  plane  for  the  vector  of  unknown 
constants  allows  the  displacement  potentials  of  the  system  to  be  evaluated  using  equations  (26), 
(55),  and  (59).  Equation  (70)  will  be  evaluated  for  only  one  forcing  function  at  a  time.  The 
coefficients  will  be  normalized  by  the  forcing  function  stress,  which  will  cast  the  output  in  terms 
of  a  transfer  function  in  the  particular  displacement,  velocity,  stress,  or  strain  of  the  cylinder. 


'-Po 

B? 

cV 

-p. 

M 

0 

(70) 


From  boundary  condition  1, 


d 


X,n^ 


SR(1,1)  -  (Xj +2pj)— 7^(p,fl)  +  -^— y^(Pja)--^/^(Pja)-A:2)ij/^(Pja),  (71) 


Si?(l,2)  =  -2lti[J/„(^ia)-^|:/„(^ia)), 


(72) 


Si?  (1,3) 


(13) 


and 


SR  (1,4)  = 


(74) 


From  boundary  condition  2, 


Si?  (2,  1) 


(75) 
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SR  {1,2)  = 


SR{2,A)  =  0. 


From  boundary  condition  3, 


5i?(3,l)  =  -/>ia)-W(Pia)  , 


5i?(3,2)  =  + 


SR  (3,  3)  =  rtH,f4-',,i('?i«)--',+  i(‘i'l‘')('^ 


SR(3,4)  =  0. 


From  boundary  condition  4, 


5i?(4,  1)  = 


SR  {A,  2)  =  -Jiq.a), 


Si?  (4,  3)  =  ikJ^^^iqyd). 


Si?  (4,  4)  = 


3  „(1)  /  N 
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DOUBLE  SOLID  CYLINDER/FLUID  INTERFACES 
BOUNDARY  CONDITIONS 

In  this  section,  we  will  formulate  the  boundary  conditions  for  the  two-cylinder  problem 
shown  in  figure  1.  At  radius  r  =  b,  the  second  cylinder/fluid  interface,  the  radial  stress  in 
cyhnder  2  is  set  equal  to  the  negative  of  the  magnitude  of  the  pressure  in  the  fluid.  Also  at  radius 
r  =  b,  the  radial  displacement  of  the  cylinder  equals  the  radial  displacement  of  the  fluid.  At 
radius  r  =  a,  there  is  continuity  of  displacements  and  stress  components  between  each  cylinder  in 
the  plane  of  contact  perpendicular  to  the  radial  coordinate,  x^^,  .  For  the  composite 

problem,  the  following  description  at  the  boundaries  exists: 


Boundary  condition  1 

C2 

=  (-Po-Ps)\ 

r  =  b  Ir  = 

(87) 

Boundary  condition  2 

C2 

II 

II 

\r  =  b^ 

(88) 

Boundary  condition  3 

C2 

^,-e 

II 

II 

<x> 

\r  =  b' 

(89) 

Boundary  condition  4 

C2 

=  b 

(90) 

r  =  b  ^  Ir 

Boundary  condition  5 

C2 

Cl 

(91) 

r  =  a 

r  =  a 

Boundary  condition  6 

C2 

Cl 

(92) 

r  -  a 

r  =  a 

Boundary  condition  7 

C2 

Cl 

(93) 

r  =  a 

r  =  a 

Boundary  condition  8 

C2 

Cl 

= 

(94) 

r  =  a 

r  =  a 

Boundary  condition  9 

C2 

Cl 

(95) 
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Boundary  condition  10 


C2 


=  w. 


Cl 


r  =  a 


(96) 


The  10  boundary  conditions  are  represented  diagrammatically  in  figure  5  for  the  double  solid 
cylinder  immersed  in  an  outer  fluid.  All  boundary  conditions  are  in  effect  over  the  entire  surface 
of  contact  at  radii  r  =  a  and  r  =  b. 


Figure  5.  Boundary  Conditions  for  the  Double  Solid  Cylinder  and  Fluid 
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SYSTEM  MATRIX 

The  previous  10  boundary  conditions  are  assembled  into  matrix  form  in  the  following 
expression: 

5(1,1)  5(1,2)  5(1,3)  5(1,4)  5(1,5)  5(1,6)  0  0  0  5(1,10)  4^^ 

5(2,1)  5(2,2)  5(2,3)  5(2,4)  5(2,5)  5(2,6)  0  0  0  0 

5(3,1)  5(3,2)  5(3,3)  5(3,4)  5(3,5)  5(3,6)  0  0  0  0  5^^ 

5(4,1)  5(4,2)  5(4,3)  5(4,4)  5(4,5)  5(4,6)  0  0  0  5(4,10) 

5(5,1)  5(5,2)  5(5,3)  5(5,4)  5(5,5)  5(5,6)  5(5,7)  5(5,8)  5(5,9)  0  cf^ 

5(6,1)  5(6,2)  5(6,3)  5(6,4)  5(6,5)  5(6,6)  5(6,7)  5(6,8)  5(6,9)  0 

5(7,1)  5(7,2)  5(7,3)  5(7,4)  5(7,5)  5(7,6)  5(7,7)  5(7,8)  5(7,9)  0 

5(8,1)  5(8,2)  5(8,3)  5(8,4)  5(8,5)  5(8,6)  5(8,7)  5(8,8)  5(8,9)  0  sf' 

5(9,1)  5(9,2)  5(9,3)  5(9,4)  5(9,5)  5(9,6)  5(9,7)  5(9,8)  5(9,9)  0  cf' 

5(10,1)  5(10,2)  0  0  5(10,5)  5(10,6)  5(10,7)  0  5(10,9)  0  M 

The  elements  of  equation  (97)  are  given  by  equations  (98)  through  (197).  From  boundary 


condition  1,  the  first  row  yields 

S(l,l)  =  +  +  +  (98) 

5(1,2)  =  +  +  +  (99) 

5  ( 1, 3)  =  -  f-|j,(<,2«) ,  (100) 

5(1,4)=  -2n,i(  r„(92«.)  -  b^yj^q^b)] ,  (101) 

j(^2^)  >  (102) 

S(l,6)  =  (103) 
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5(1,7)  =  0, 

5(1,8)  =  0, 

5(1,9)  =  0, 

5(1,10)  =  co2p^M//y^g2^)  . 


(104) 

(105) 

(106) 

(107) 
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From  boundary  condition  2, 


5(2,  1)  =  (108) 

5(2,2)  =  /2ii2^|-7>2^),  (109) 

iVi^kn 

5(2,3)  =  (110) 

5(2,4)  =  -^Y^iq^b),  (111) 


5(2,5)  (1,2) 

5(2,6)  =n^(r„^,(,^i,)(!!^-i2^-|.y„,,(„,i,)(l^)-^y,^, (,,,(,)),  (ii3) 


5(2,7)  =0,  (114) 

5(2,8)  =0,  (115) 

5(2,9)  =0,  (116) 

and 

5(2, 10)  =0.  (117) 
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From  boundary  condition  3, 

S(3,l)  = 

5(3,2)  =  <“’* 

5(3,3)  (120) 

5(3,4)  =|42(-^I'„(42«’)  +  ^|:W)-^2W)).  (121) 

5(3,5)  =  i<:H2(^-(„*l(?2'’)--(»+l(«2(’>(^))'  (122) 

5(3,6)  =  ;<rH2(|:l'„„(?2fc)-l'„  +  i(92(’)(^)).  (123) 

S(3,7)  =0,  (124) 

5(3,8)  =0,  (125) 

5(3,9)  =0,  (126) 

and 

5(3, 10)  =0.  (127) 


30 


TR  11,043 


From  boundary  condition  4, 


5(4,1) 

=  |:W). 

(128) 

5(4,  2) 

(129) 

5(4,3) 

(130) 

5(4,4) 

(131) 

5(4,5) 

=  ikJ^^^iq^b), 

(132) 

5(4,  6) 

(133) 

5(4,7) 

=  0, 

(134) 

5(4,8) 

=  0, 

(135) 

5(4,9)  =  0, 
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From  boundary  condition  5, 

S(5,l)  =  +  +  + 

S(5,2)  =  (?l2  +  2^l2)^l^>2«)  +  7|:^«(^2«)-V>2^)(^  +  ^^)’ 
S  (5,  3)  =  -2|X2^(-^„(^2^)  -  ’ 

S  (5,  4)  =  -2H2^(  ’ 

5(5,5)  =  i2\i2k^J„^M^), 

5(5,6)  =  /2^2^^r„^i(^2«)’ 

5(5,7)  =  (?li+2Hi)^/„(Pia)  +  ^|:4(Pi«)-V>i«)(^  +  ^^)’ 
5(5,8)  =  — 

5(5,9)  =  -2Fii^^4+i(^i'')’ 

and 

5(5, 10)  =  0. 


(138) 

(139) 

(140) 

(141) 

(142) 

(143) 

(144) 

(145) 

(146) 

(147) 
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From  boundary  condition  6, 


5(6,1)  =  i2\i2k-^J„(P2a), 


5(6,2)  = 


i[ir,kn 

5(6,3)  =  -^JJiq^a), 


iii^kn 

5(6,4)  =  -^Y^iq^a), 


(148) 

(149) 

(150) 

(151) 


5(6,5)  -  +  n  J  j  ’ 


V 


a  )  9r2 


5(6,6)  -  ^^n  +  l(‘3'2^\  ^  J 


5(6,7)  =  -il]i^k^JJ^^a), 
i\x.,kn 

5(6,8)  = - ^7„(<7i«), 


n  +  \ 


5(6,9)  -  l(^l^)(^  ^2  ^  + 


n  +  1  ^  d' 


(154) 


(155) 


(156) 


and 


5(6, 10)  =  0. 


(157) 
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From  boundary  condition  7, 

s  (7, 2)  =  -^l^iYJp,a)-^Y„(p,a) J ,  (159) 

5(7,3)  +  (160) 

+  (161) 

5(7,5)  =  rtn,(l7„,,(?j«)-/„,,(«2«)(i^)),  (162) 

5(7,6)  =  »H2(|:1'„,,(«2«)-1'„,(92«)(^)).  (163) 

5(7,7)  =-?^(/^(;,,<,)-a|-/>,a)),  (164) 

5(7,8)  =-n,(-^/„(„,a)  +  i|7^(,,a)-i^7„(,,a)),  (165) 

5(7,9)  =-rtH,(|:7,.,(9,<i)-7„,(«iO)(i^)).  (166) 

and 

5(7,  10)  =0.  (167) 
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From  boundary  condition  8, 


(168) 

5(8,2)  =  yj.p^a). 

(169) 

5(8,3)  = 

(170) 

5(8,4)  =  -F„(^2«), 

(171) 

5(8,5)  =  ikJ^^^iq^a), 

(172) 

5(8,6)  =  ikY^^^{q^a), 

(173) 

5(8,7)  =-|/>,a). 

(174) 

5(8,8)  =-^7„(?,a). 

(175) 

5  (8,  9)  =  j(gja), 

(176) 

5(8,  10)  =  0. 

(177) 
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From  boundary  condition  9, 


5(9,1) 

(178) 

5(9,  2) 

(179) 

5(9,3) 

(180) 

5(9,4) 

(181) 

5(9,5) 

(182) 

5(9,  6) 

=  ikY^^^iq^a), 

(183) 

5(9,  7) 

(184) 

5(9,  8) 

(185) 

5(9,  9) 

=  -ikJ^^^iq^a), 

(186) 

and 

S{9,  10)  =0.  (187) 
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From  boundary  condition  10, 


5(10,1)  =  ikJ^ip^a),  (188) 

5(10,2)  =  ikYJ^^a),  (189) 

5(10,3)  =0,  (190) 

5(10,4)  =0,  (191) 

5(10,5)  ^•^n+ ’  (192) 

S(io,6)  («,„),  (193) 

5(10,7)  ^  -ikJ^ip^a),  (194) 

5(10,8)  =0,  (195) 

5(10,9)  =  ^•^„  +  i(<?i«)  +  ^^-^/„  +  i(9ia),  (196) 

and 

5(10,10)  =0.  (197) 
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RESULTS 

BARE  OPTICAL  FIBER 
Optical  Phase  Sensitivity 

The  first  numerical  example  will  be  a  calculation  using  the  solid  cylinder  model  with 
properties  of  a  bare  optical  glass  fiber  immersed  in  water  (tables  1  and  2).  The  optical  phase 
sensitivity  ^  will  be  computed  using  the  following  expression,^**  where  p  =  P^,  or  Pq: 

2 

^  .  (198) 


Table  1.  Optical  Fiber  Properties 


Property 

Definition 

E.  =  7.2xl0‘°^ 

Young’s  Modulus 

=  0.005 

Structural  Loss  Factor 

p,  =  2600li 

HI'’ 

Density 

Vj  =  0.17 

Poisson’s  Ratio 

a  =  62.5xlO"V 

Cladding  Outer  Radius 

=  4.5x10  ^m 

Core  Radius 

rj  =  1.0x10“% 

Calculation  Radius 

=  1.46 

Refractive  Index 

Pjj  =  0.126 

Pockel  Coefficient 

Pj2  =  0-27 

Pockel  Coefficient 
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Table  2.  Outer  Fluid  Properties 


Property 

Definition 

=  1000.0^ 

Density 

c  =  1500.0— 
sec 

Speed  of  Sound 

As  was  previously  discussed,  the  system  matrix,  equation  (70),  is  evaluated  for  the  optical 
fiber  properties  listed  above.  The  calculation  radius  is  near  to  the  center  of  the  optical  fiber  core 
at  Tj  =  1.0xl0"^m.  The  coefficients  obtained  from  the  matrix  inversion  normalized  by  the 
forcing  function  stress  are  substituted  into  equation  (26)  for  the  displacement  potentials,  and  next 
into  equations  (43)  and  (45)  for  the  strains.  The  strains  are  then  used  to  obtain  the  optical  phase 
sensitivity  via  equation  (198). 


Numerical  Considerations 

The  procedure  just  described  is  repeated  over  the  domain  of  wavenumber  and  frequency  at 
regular  intervals.  For  the  «  =  0  case,  the  frequency  range  is  1.0  to  500.0  Hz  at  intervals  of  2.5 
Hz,  and  the  wavenumber  range  is  -2.0  to  2.0  rad/m  at  intervals  of  0.02  rad/m.  These  ranges 
provide  an  image  size  of  200-by-200  points  before  interpolation  for  display.  For  the  n  —  1  case, 
the  frequency  range  remains  the  same;  however,  the  wavenumber  range  is  increased  to  -200.0  to 
200.0  rad/m.  The  wavenumber  interval  is  increased  so  that  the  final  calculated  image  size  is 
consistent  with  that  of  the  n  =  0  case.  A  numerical  artifact  is  apparent  in  the  simulations  due  to 
the  finite  wavenumber  and  frequency  interval  when  sharp  resonances  are  present.  This  is  most 
evident  in  the  n  =  0  simulations  described  later  in  the  report.  The  extensional  wave  resonance  is 
not  continuously  resolved,  but  appears  as  a  series  of  dots  in  both  the  magnitude  and  phase 
surfaces.  Small  asymmetries  with  respect  to  wavenumber  are  also  evident  in  the  magnitude  cuts 
as  well.  The  first  problem  is  a  function  of  the  graphics  software  used  to  display  the  wire  frame 
surfaces,  and  the  second  problem  results  when,  for  a  given  frequency,  slightly  different 
wavenumber  locations  are  being  calculated  in  the  minus  and  plus  wavenumber  half¬ 
planes.  These  problems  could  be  eliminated  by  decreasing  both  the  frequency  and  wavenumber 
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interval  by  several  orders  of  magnitude.  A  corresponding  increase  in  calculation  time  would 
result.  The  200-by-200  point  image  (40,000  points)  was  selected  as  an  appropriate  compromise 
between  calculation  time  and  the  requirement  to  convey  the  physics  embodied  in  the  simulations. 

Evaluation  of  the  equations  required  the  use  of  the  quad  precision  complex  data  type.  The 
language  of  choice  is  FORTRAN  (version  3.0)  running  on  a  Sun  Microsystems  SPARC  Station 
model  lOSX-512,  Operating  System  Solaris  2.3.  In  the  course  of  this  research,  the  author 
discovered  a  “bug”  in  the  Fortran  complier  when  using  quad  precision  complex  functions. 
Therefore,  calculations  were  limited  to  double  precision  complex.  However,  even  with  the 
accuracy  afforded  with  the  double  precision  data  type,  numerical  accuracy  problems  are 
encountered  with  the  bare  fiber  parameters  and  the  n  =  1  excitation.  As  shown  later  in  the  report, 
numerical  accuracy  problems  will  manifest  themselves  as  random  spots  concentrated  at  low 
frequency  (below  50  Hz)  for  all  wavenumbers  and  at  all  frequencies  for  wavenumbers  below 
approximately  +50  rad/m.  In  this  case,  the  occurrence  of  such  problems  is  due  to  the  choice  of 
optical  fiber  material  properties  and  dimensions.  Several  other  cases  (not  shown)  were 
evaluated,  one  being  a  urethane  rod  (approximating  a  solid  garden  hose),  which  did  not  display 
any  of  the  numerical  accuracy  problems  evident  in  the  optical  fiber  simulations  presented. 

The  FORTRAN  algorithms  developed  to  produce  the  simulations  contained  here  are 
published  in  a  separate  document,  which  will  be  available  (for  one  year  only  from  the  date  of 
publication  of  this  technical  report)  via  e-mail  at  peloquin@cascade.nl.nuwc.navy.mil. 


CIRCUMFERENTIAL  ORDER  SHAPES 

For  the  simulations  that  follow  in  the  remainder  of  this  report,  the  motion  of  the  cylinder  will 
take  the  form  depicted  in  one  of  the  two  circumferential  order  shapes  shown  in  figures  6  and  7. 
Figure  6  depicts  the  longitudinal  x-displacement  that  occurs  in  the  cylinder  as  energy  is 
propagating  in  the  first  branch  of  wave  propagation  corresponding  to  circumferential  order 
number  n  =  0.  This  first  branch  of  wave  propagation  is  commonly  referred  to  as  an  extensional 
wave. 
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The  displacement  shape  depicted  in  figure  7  causes  the  cylinder  to  undergo  bending  in  the 
rx-plane.  The  nonaxisymmetric  excitations  considered  here  excite  this  mode  of  wave 

propagation.  As  energy  is  propagating  in  the  branches  of  the  n  =  1  simulations  shown,  the 
cylinder  is  undergoing  a  vibration  principally  depicted  by  the  deformation  shape  of  figure  7. 

The  cylinder  displacements  in  figures  6  through  8  have  been  magnified  for  the  purpose  of 
visualizing  the  deformation.  It  should  be  remembered  that  this  is  a  small  strain  analysis,  suitable 
for  analyzing  problems  undergoing  loading  on  an  acoustic  scale. 
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Figure  6.  Circumferential  Order  Number  n  =  0,  First  Branch 


Figure  7.  Circumferential  Order  Number  n  =  1,  First  Branch 
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Transfer  Surfaces 

The  simulation  for  axisymmetric  radial  pressure  excitation  results  in  transfer  surfaces  for  each 
strain  component  and  the  optical  phase  sensitivity,  as  displayed  in  figures  9  through  16.  The  cuts 
through  the  magnitude  surfaces  show  the  character  of  the  field  quantities  as  functions  of 
wavenumber.  The  wavenumber  axis  is  marked  with  major  tick  marks  at  the  limits  of  the  range 
for  acoustic  plane  wave  propagation  in  water.  These  points  correspond  to  minus  and  plus 
endfire.  The  predominant  motion  of  the  cylinder  is  in  the  longitudinal  x-direction,  as  depicted  in 
figure  6.  The  major  feature  in  this  simulation  is  the  supersonic  extensional  wave  at  557 1  m/sec. 


Since  the  extensional  wave  phase  velocity  is  supersonic,  it  generates  a  radial  component  of 
cylinder  displacement  that  can  be  projected  onto  the  fluid  propagation  direction;  energy  will  thus 
be  radiated  into  the  outer  fluid  at  the  fluid  dilatational  wave  speed.  Note  that  a  projection  of  the 
extensional  wavelength  of  the  cylinder,  L^,  onto  the  corresponding  fluid  wavelength,  L„  exists. 
Constant  pressure  phase  fronts  in  the  fluid  will  leave  the  cylinder  at  oc  =  74  degrees  relative  to  the 


Figure  8.  Radiation  of  Pressure  From  the  Cylinder  Surface  at  105  Degrees 
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longitudinal  x-axis  (figure  8),  the  fluid  medium  being  water  with  =  1500  m/sec,  as  given  in  the 
list  of  outer  fluid  properties.  Figure  8  represents  the  qualitative  relationship  that  exists  between 
the  wavelengths  for  every  location  in  the  wavenumber-frequency  plane  along  the  branch 
representing  extensional  wave  propagation. 

Conversely,  a  pressure  wave  in  the  outer  fluid  incident  upon  the  cylinder  at  74  degrees,  as 
shown  in  figure  8,  would  excite  an  extensional  wave  response  in  the  cylinder.  Remembering 
that  kg  =  InlLg  and  kg  =  co/c^,  the  relationship  between  the  incident/radiated  angle,  wavelength, 
phase  velocity,  and  wavenumber  is 


L  c  k 

a  =  acos —  =  acos—  =  acos-^. 


(199) 


To  review,  in  an  infinite  isotropic  elastic  media,  the  dilatational  wave  phase  velocity  is  given 
by  equation  (200)  as 


fp  11-hvJ(1-2vJ 


(200) 


For  a  dilatational  wave,  the  particle  velocity  and  the  wave  propagation  vector  are 
coincident.  For  the  wavenumber  and  frequency  range  shown  in  the  simulations,  the 
manifestation  of  this  wave  in  the  cylinder  produces  a  modified  phase  velocity  known  as  the  bar 
velocity  or  extensional  wave  phase  velocity,  given  by  equation  (201)  as 


c 


e 


(201) 


The  simulation  for  axisymmetric  longitudinal  shear  stress  results  in  figures  17  through  24. 
Note  again  the  prominent  extensional  wave  resonant  peak.  The  strain  and  optical  phase 
sensitivity  (in  the  fiber)  to  longitudinal  shear  stress  is  some  30  to  110  dB  greater  than  it  is  for 
sensitivity  to  radial  pressure  excitation.  This  is  evident  by  comparing  the  difference  in  levels 
between  figures  10  and  18,  figures  12  and  20,  figures  14  and  22,  and  figures  16  and  24.  At  the 
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extensional  wave  resonance  peak,  the  excess  is  110  dB.  At  the  edge  of  the  acoustic  cone,  the 
excess  is  70  dB,  and  at  broadside  {k  =  0),  the  excess  is  30  dB.  Thus,  it  is  quite  clear  that 
sensitivity  to  longitudinal  shear  stress  exceeds  sensitivity  to  radial  pressure  excitation  by  over  5 
orders  of  magnitude. 

Two  nonaxisymmetric  simulations  have  been  performed  for  circumferential  order  number 
n=  Pg,  and  The  nonaxisymmetric  radial  pressure  simulation  is  contained  in  figures  25 
through  32.  This  excitation  imparts  a  bending  excitation  to  the  fiber  whose  motion  exists 
primarily  in  the  rx-plane.  This  excitation  results  in  dispersive  wave  propagation,  which 
produces  the  nonstraight-line  characteristic  of  the  branches  in  figure  25.  Energy  propagates  at  a 
much  slower  rate  in  the  n  =  1  branches  than  in  the  n  =  0  branches  seen  in  the  earlier  figures.  The 
phase  velocity  ranges  from  12  to  20  m/sec  over  the  frequency  range  from  1  to  500  Hz,  as  shown 
in  figure  25. 

An  R  =  1  longitudinal  shear  stress  excitation  produces  a  deformation,  principally  in  the 
rx-plane  of  the  cylinder.  The  strains  and  optical  phase  sensitivities  for  this  simulation  are  shown 
in  figures  33  through  40.  For  the  n  =  1  case,  the  magnitude  of  the  strains  produced  by  the 
longitudinal  shear  stress  excitation  is  20  to  100  dB  smaller  than  the  magnitude  of  the  strains 
produced  by  the  radial  pressure  excitation  (compare  figure  34  with  26).  Both  of  the  n  =  \ 
simulations  suffer  from  the  numerical  accuracy  problem  described  earlier,  which  is  the  reason  for 
the  spottiness  in  the  r  =  1  figures.  All  of  the  r  =  1  simulations  are  performed  at  a  calculation 
radius  of  3.936  x  10~^  meters,  which  is  outside  of  the  fiber  core  but  within  the  fiber  outer 
diameter. 

Mechanical  damping  in  the  cylinder  is  incorporated  by  the  use  of  the  structural  loss  factor. 
This  structural  loss  factor  is  the  imaginary  component  of  the  complex  modulus  of  elasticity: 


E*j  =  E.{\+i(,p  ,  (202) 

where  j  equals  1  or  2,  corresponding  to  the  first  and  second  layer,  respectively.  It  is  therefore 
possible,  within  the  formulation  presented  here,  to  have  different  structural  loss  factors  for  each 
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cylinder.  This  is,  in  fact,  how  the  double  layer  simulations  were  performed. 

When  phase  is  used  as  a  label  for  the  lower  of  the  two  images  in  the  figures  that  follow,  it  is 
the  phase  between  the  applied  excitation  stress  and  the  particular  variable  under  consideration. 
This  is  the  phase  between  the  applied  stress  and  response  of  the  cylinder.  It  should  not  be 
confused  with  the  optical  phase  sensitivity  defined  by  equation  (198). 
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Axisymmetric  Radial  Pressure  Excitation  (n  —  0) 
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Figure  9.  Radial  Strain  Transfer  Surface:  Magnitude  =  lOLog(e^Po)^  and  Phase  in 

Degrees  (One  Layer, «  =  0) 


Figure  10.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  9 
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Figure  13.  Longitudinal  Strain  Transfer  Surface:  Magnitude  =  IQiLo^iz^jP and  Phase 

in  Degrees  (One  Layer, «  =  0) 
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Figure  14.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  13 
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Figure  15.  Optical  Phase  Sensitivity  Strain  Transfer  Surface:  Magnitude  =  10Log(A(|)/(])/P^) 

and  Phase  in  Degrees  (One  Layer,  w  =  0) 


Figure  16.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  15 
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Axisymmetric  Longitudinal  Shear  Stress  Excitation  (n  —  0) 
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Figure  17.  Radial  Strain  Transfer  Surface:  Magnitude  =  10Log(ejPx)^  and  Phase  in 

Degrees  (One  Layer,  n  =  0) 


Figure  18.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  17 
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Figure  19.  Circumferential  Strain  Transfer  Surface:  Magnitude  =  lOLog(800/Pj.)^  and 

Phase  in  Degrees  (One  Layer,  n  =  0) 
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Figure  20.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  19 
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Figure  21.  Longitudinal  Strain  Transfer  Surface:  Magnitude  =  10Log(e;,.;^/P;^)^  and  Phase 

in  Degrees  (One  Layer,  n  =  0) 


Figure  22.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  21 
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Figure  23,  Optical  Phase  Sensitivity  Strain  Transfer  Surface:  Magnitude  =  10Log(A(j)/(j)/PD 

and  Phase  in  Degrees  (One  Layer,  n  =  0) 
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Figure  24.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  23 


Figure  25.  Radial  Strain  Transfer  Surface:  Magnitude  =  10Log(£r/Po)^  and  Phase  in 

Degrees  (One  Layer,  u  =  1) 


Figure  26.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  25 
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27.  Circumferential  Strain  Transfer  Surface:  Magnitude  =  lOLog(800//*^)^  and 
Phase  in  Degrees  (One  Layer,  n  =  1) 
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Figure  28.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  27 
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Figure  29.  Longitudinal  Strain  Transfer  Surface:  Magnitude  =  and  Phase 

in  Degrees  (One  Layer,  n  =  1) 
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Figure  30.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  29 
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.  Optical  Phase  Sensitivity  Strain  Transfer  Surface:  Magnitude  =  10Log(A(|)/(|)/P^) 
and  Phase  in  Degrees  (One  Layer,  /i  =  1) 


Figure  32.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  31 
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Nonaxisymmetric  Longitudinal  Shear  Stress  Excitation  (n  =  1) 
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Figure  33.  Radial  Strain  Transfer  Surface:  Magnitude  =  10Log(8;./F;,)^  and  Phase  in 

Degrees  (One  Layer,  n  =  1) 


Figure  34.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  33 
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Figure  35.  Circumferential  Strain  Transfer  Surface:  Magnitude  =  lOLog(8e0/P^)^  and 

Phase  in  Degrees  (One  Layer, «  =  1) 
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Figure  36.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  35 
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17.  Longitudinal  Strain  Transfer  Surface:  Magnitude  =  10Log(e;ex/^jc)^  and  Phase 

in  Degrees  (One  Layer,  n  =  1) 
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Figure  38.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  37 
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Effects  of  Poisson  Ratio 

The  Poisson  ratio  influences  the  shape  of  the  dynamic  response  as  a  function  of  wavenumber. 
Figure  41  compares  10Log(8^,/Po)^  for  three  different  values  of  Poisson  ratio.  The  case  of  =  0 
results  in  the  flat  response  indicated  by  the  solid  line;  the  extensional  wave  is  not  excited  for  this 
case.  When  =  0,  the  three  strains  are  decoupled  from  each  other.  Since  the  excitation,  P is 
in  the  radial  direction  there  is  no  way  to  perturb  the  longitudinal  direction,  hence  the  lack  of 
extensional  wave  response.  The  Poisson  ratio  increasing  from  0.17  to  0.45  results  in  an 
increasingly  large  extensional  wave  resonance  amplitude.  As  the  value  of  the  Poisson  ratio 
increases,  the  relative  difference  between  the  /:  =  0  value  of  10Log(e^/Po)^  fo®  value  of 
lOLog(e^^P^)^  above  the  extensional  wave  resonance  (flat  plateau)  increases.  This  increase  of 
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the  difference  in  relative  levels  above  and  below  extensional  wave  resonance  is  also  evident  when 
a  coating  is  applied  to  the  fiber,  which  will  be  seen  in  the  next  section. 

COATED  OPTICAL  FIBER 

The  bare  optical  fiber  of  the  previous  section  is  coated  with  an  ALCRYN  layer.  ALCRYN  is  a 
viscoelastic  rubber  whose  properties  are  listed  in  table  3.  The  dynamic  material  properties  of 


Table  3.  Optical  Fiber  and  ALCRYN  Properties 


Optical  Fiber  Property 

ALCRYN  Property 

Definition 

E,  =  7.2x10*°-^ 

7  N 

E.  =  1.64x10  ^ 

Young’s  Modulus 

I! 

O 

o 

o 

Cj  =  0.47 

Structural  Loss  Factor 

p,  =  2600^ 

P2  =  1000^ 

Density 

Vj  =  0.17 

Vj  =  0.45 

Poisson’s  Ratio 

a  =  62.5xl0~^m 

Cladding  Outer  Radius 

=  4.5x10  ^m 

Core  Radius 

Tj  =  1.0xl0~^m 

Calculation  Radius 

a  =  62.5x10  ^m 

Coating  Inner  Radius 

b  =  4.0xl0~^m 

Coating  Outer  Radius 

=  1.46 

Refractive  Index 

Pjj  =  0.126 

Pockel  Coefficient 

=  0.27 

Pockel  Coefficient 
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ALCRYN  were  measured  with  a  Metravib  dynamic  mechanical  analyzer  at  the  Naval  Undersea 
Warfare  Center  Detachment  in  New  London,  Connecticut.  The  cladding  outer  radius  is  always 
equal  to  the  inner  radius  of  the  coating  (see  table  3),  consistent  with  boundary  condition  8. 


Transfer  Surfaces,  Two  Layers 

Similar  to  the  previous  section  for  the  bare  optical  fiber,  the  response  of  the  coated  fiber  will 
be  simulated  using  the  two-layer  infinite  cylinder  model  for  two  different  circumferential  order 
numbers  and  three  different  excitation  types. 

The  response  to  the  first  excitation,  is  covered  in  figures  42  through  49.  Each  strain 
component  is  displayed  individually  along  with  the  corresponding  optical  phase  sensitivity  image 
(figure  48)  and  the  corresponding  cuts  through  the  image  (figure  49).  The  most  dramatic  effect 
of  the  coating  is  to  slow  the  extensional  wave  to  190  m/sec,  down  from  5,571  m/sec  for  the  bare 
fiber  of  figure  9.  In  the  radial  strain  simulation  at  40  Hz  (figure  43),  the  extensional  wave 
resonance  peak  is  greater  than  60  dB  above  the  k  =  0  level.  The  difference  grows  to  over  80  dB 
in  the  optical  phase  sensitivity  cuts  shown  in  figure  49.  The  effects  of  material  properties  (first 
observed  in  the  case  of  the  bare  fiber),  which  create  a  relative  difference  in  strain  response 
between  k  =  0  and  the  wavenumber  region  at  and  above  extensional  wave  resonance,  are  much 
greater  for  the  coated  fiber  case.  Additionally,  since  the  coating  is  much  softer  than  the  fiber,  it 
contributes  mass  to  the  composite,  adding  very  little  stiffness.  This  is  why  the  extensional  wave 
phase  velocity  slowed  so  dramatically. 

The  simulations  for  the  response  to  longitudinal  shear  stress  excitation  are  shown  in  figures  50 
through  57.  Similar  to  the  excitation  (figure  42),  the  extensional  wave  phase  velocity  has 
slowed  to  190  m/sec.  Comparison  is  made  between  the  optical  phase  response  to 
figure  58.  Notice  the  tremendous  increase  in  sensitivity  to  shear  stress  P^  over  the  sensitivity  to 
normal  pressure  P^ — on  the  order  of  50  to  70  dB. 

All  of  the  nonaxisymmetric  excitations,  P^,  P^,  and  Pq,  produce  dispersive  branches  of  wave 
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propagation,  as  seen  in  figures  59  through  82.  The  motion  of  the  cylinder  follows  the 
deformation  shape  depicted  in  figure  7.  Excitations  and  produce  strains  of  approximately 
equivalent  magnitude.  Excitation  produces  strains  some  10  dB  lower  than  and  Pq,  as  can 
be  observed  by  inspection  of  these  figures. 
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Axisymmetric  Radial  Pressure  Excitation  (n  =  0),  Two  Layers 
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Figure  42.  Radial  Strain  Transfer  Surface:  Magnitude  =  10Log(e^/Po)^  Phase  in 

Degrees  (Two  Layers,  n  =  0) 
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Figure  43.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  42 
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Figure  44.  Circumferential  Strain  Transfer  Surface:  Magnitude  =  lOLogCsee/P^)^  and 

Phase  in  Degrees  (Two  Layers,  n  =  0) 


Figure  45.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  44 
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Figure  46.  Longitudinal  Strain  Transfer  Surface:  Magnitude  =  \QiLo%izjPof  and  Phase  in 

Degrees  (Two  Layers,  n  =  0) 
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Figure  47.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  46 
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Figure  48.  Optical  Phase  Sensitivity  Strain  Transfer  Surface:  Magnitude  =  10Log(A(j)/(|)/Pp) 

and  Phase  in  Degrees  (Two  Layers, «  =  0) 


Figure  49.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  48 
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Axisymmetric  Longitudinal  Shear  Stress  Excitation  (n  =  0),  Two  Layers 
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Figure  50.  Radial  Strain  Transfer  Surface:  Magnitude  =  ItSLogizjPx?  and  Phase  in 

Degrees  (Two  Layers,  n  =  0) 


Figure  51.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  50 
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Figure  52.  Circumferential  Strain  Transfer  Surface:  Magnitude  =  lOLogCeee/P;^)^  and 

Phase  in  Degrees  (Two  Layers,  n  =  0) 
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Figure  53.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  52 
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Figure  54.  Longitudinal  Strain  Transfer  Surface:  Magnitude  =  Phase 

in  Degrees  (Two  Layers,  n  =  0) 
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Figure  55.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  54 
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Figure  56.  Optical  Phase  Sensitivity  Strain  TVansfer  Surface:  Magnitude  =  10Log(A(|)/(j)/Pjf) 

and  Phase  in  Degrees  (Two  Layers,  n  =  0) 
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Figure  57.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  56 
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Nonaxisymmetric  Radial  Pressure  Excitation  (n  =  l),  Two  Layers 
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Figure  59.  Radial  Strain  Transfer  Surface:  Magnitude  =  10Log(8r/Po)^  and  Phase  in 

Degrees  (Two  Layers,  n  =  1) 
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Figure  60.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  59 
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61.  Circumferential  Strain  Transfer  Surface:  Magnitude  =  lOLog(ee0/F<,)^  and 
Phase  in  Degrees  (Two  Layers,  n  =  1) 
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Figure  62.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  61 
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Figure  63.  Longitudinal  Strain  IVansfer  Surface:  Magnitude  =  lOLogCe^^P^)^  and  Phase  in 

Degrees  (Two  Layers,  n  =  1) 


MAGNITUDE 


0 

WAVENUMBER  (rad/m) 


200 


Figure  64.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  63 
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»5.  Optical  Phase  Sensitivity  Transfer  Surface:  Magnitude  =  10Log(A(t)/(j)/P and 
Phase  in  Degrees  (Two  Layers,  n  =  1) 
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Figure  66.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  65 
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Nonaxisymmetric  Longitudinal  Shear  Stress  Excitation  (n  =  1),  Two  Layers 


Figure  67.  Radial  Strain  Transfer  Surface:  Magnitude  =  10Log(e;.;yP;^)^  and  Phase  in 

Degrees  (Two  Layers,  n  =  1) 


Figure  68.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  67 


81 


lOLog(800/P^)^  j,  WAVENUMBER  (rad/m)  WAVENUMBER  (rad/ 


69.  Circumferential  Strain  Transfer  Surface:  Magnitude  =  lOLog(80e/P^)^  and 
Phase  in  Degrees  (Two  Layers, «  =  1) 


Figure  70.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  69 
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Figure  71.  Longitudinal  Strain  Transfer  Surface:  Magnitude  =  10Log(e_^//*_j.)^  and  Phase  in 

Degrees  (Two  Layers,  n  =  1) 


Figure  72.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  71 
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Figure  73.  Optical  Phase  Sensitivity  Strain  Transfer  Surface:  Magnitude  =  10Log(A<t)/(|)/P^)^ 

and  Phase  in  Degrees  (Two  Layers,  n  =  1) 
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Figure  74.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  73 
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Nonaxisymmetric  Circumferential  Shear  Stress  Excitation  (n- 1),  Two  Layers 


Figure  75.  Radial  Strain  IVansfer  Surface:  Magnitude  =  lOLog(e;./P0)^  and  Phase  in 

Degrees  (Two  Layers,  n  =  l) 
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Figure  76.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  75 
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!  77.  Circumferential  Strain  Transfer  Surface:  Magnitude  =  10Log(eee/P e)^  and 
Phase  in  Degrees  (Two  Layers,  n  =  1) 


Figure  78.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  77 
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Figure  79.  Longitudinal  Strain  TVansfer  Surface:  Magnitude  =  lOLog(8j^//*0)^  and  Phase  in 

Degrees  (Two  Layers,  n  =  1) 


Figure  80.  Frequency  Cut  at  453.19  Hz  Through  the  Magnitude  Surface  of  Figure  79 
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EFFECTS  OF  THE  SECOND  LAYER 

We  will  now  consider  the  effects  of  a  coating  on  the  optical  phase  sensitivity  of  a  bare  and 
coated  optical  fiber  subject  to  a  harmonic  radial  pressure  excitation  Pg.  The  application  of  a 
coating  results  in  an  effect  that  is  similar  to  what  occurs  when  there  is  a  variation  in  Poisson  ratio 
(figure  41). 

Figure  83  compares  the  optical  phase  sensitivity  for  a  bare  and  coated  optical  fiber  by 
overlaying  cuts  through  figures  15  and  48  at  40  Hz.  The  wavenumbers  corresponding  to 
acoustic  propagation  speeds  for  directions  parallel  to  the  longitudinal  axis  of  the  fiber  are 
indicated  on  figure  83  for  both  air  and  water.  The  foremost  effect  of  the  coating  is  to  increase  the 
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Optical  phase  sensitivity  of  the  optical  fiber  in  the  region  of  the  extensional  wave  resonance.  The 
optical  phase  sensitivity  is  much  lower,  approximately  80  dB,  in  the  wavenumber  region 
corresponding  to  that  of  the  bare  optical  fiber  extensional  wave.  The  choice  of  optical  phase 
sensitivity  for  this  discussion  is  arbitrary;  the  comments  are  equally  valid  for  any  one  of  the  strain 
components  since  the  optical  phase  sensitivity  depends  on  the  strains  as  given  by  equation  (198). 
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CONCLUSIONS 

From  first  principles,  in  the  field  of  dynamic  elasticity,  a  model  has  been  derived  for  the  steady 
state  response  of  a  two-layer  infinite  elastic  solid  cylinder  in  contact  with  an  infinite  outer  fluid 
subject  to  forced  harmonic  excitation  at  the  solid/fluid  interface.  The  cylinders  are  treated  as 
viscoelastic  solids,  and  the  fluid  is  treated  as  an  inviscid  ideal  fluid.  Both  axisymmetric  and 
nonaxisymmetric  excitations  are  included  in  the  analysis. 

The  solution  is  a  closed  form  analytical  result  to  the  point  of  assembly  of  the  dynamic  system 
matrix.  The  system  matrix  is  evaluated  numerically  at  each  discrete  point  in  the  wavenumber- 
frequency  domain  for  the  field  parameters,  displacements,  strains,  stresses,  velocities,  and  fluid 
pressure.  The  elasticity  model  used  for  the  cylinder  is  quite  robust  and  will  accurately  describe 
the  behavior  of  a  wide  class  of  materials  whose  mechanical  properties  can  be  described  by  two  of 
the  following  macroscopic  material  properties:  Young’s  modulus,  Poisson  ratio,  or  shear 
modulus.  These  materials  range  from  closed-cell  foam  to  optical  glass  fiber. 

The  simulations  presented  are  for  the  single  cylinder  with  an  outer  fluid  and  a  double  cylinder 
with  an  outer  fluid.  The  simulations  for  the  single  cylinder  were  performed  using  material 
properties  consistent  with  those  of  optical  glass  fiber.  Two  excitations  and  were  used  for 
these  simulations  at  two  different  circumferential  order  numbers  (n  =  0  and  «  =  1).  The  first 
branch  of  wave  propagation  supported  by  the  single  cylinder  of  optical  glass  fiber  is  the 
extensional  wave,  with  energy  propagating  at  a  phase  velocity  of  5571  m/sec.  From  these 
simulations,  it  is  evident  that  the  strain  and  optical  phase  sensitivity  (in  the  fiber)  to  longitudinal 
shear  stress  is  30  to  110  dB  greater  than  the  fiber  sensitivity  to  radial  pressure  excitation.  At  the 
extensional  wave  resonance  peak,  the  shear  stress  sensitivity  is  in  excess  of  110  dB.  At  the  edge 
of  the  acoustic  cone,  the  shear  stress  sensitivity  is  in  excess  of  70  dB,  and  at  broadside  {k  =  0),  this 
excess  is  30  dB.  This  analysis  states  that  the  sensitivity  of  the  optical  fiber  to  shear  stress  is  five 
orders  of  magnitude  higher  than  its  sensitivity  to  radial  pressure.  Therefore,  the  radial  pressure 
response  will  be  dominated  by  the  response  to  shear  stress. 

Two  nonaxisymmetric  simulations  have  been  performed  for  circumferential  order  number 
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n  =  1:  Po  and  P^.  These  excitations  result  in  branches  of  dispersive  wave  propagation.  Energy 
propagating  in  the  n  =  1  branches  travels  at  a  much  slower  rate  than  energy  propagating  in  the 
n  =  0  branches.  For  the  n  =  1  case,  the  phase  velocity  ranges  from  12  to  20  m/sec  over  the 
frequency  range  from  1  to  500  Hz.  The  n  =■  \  excitations  produce  a  deformation,  principally  in 
the  rx-plane  of  the  cylinder,  whereas  the  principal  deformation  of  the  n  =  0  excitation  is  in  the 
longitudinal  jc-direction.  For  the  n  =  1  case,  the  magnitude  of  the  strains  produced  by  the 
longitudinal  shear  stress  excitation  is  20  to  100  dB  smaller  than  the  magnitude  of  the  strains 
produced  by  the  n  =  1  radial  pressure  excitation. 

The  simulations  performed  with  the  double  cylinder  model  were  used  to  analyze  the  dynamic 
response  of  the  optical  glass  fiber  coated  with  ALCRYN.  The  dynamic  material  properties  for 
the  Young’s  modulus  of  ALCRYN  were  measured  with  a  Metravib  dynamic  mechanical  analyzer 
at  the  Naval  Undersea  Warfare  Center  Detachment  in  New  London,  Connecticut.  The  and  P^ 
excitations  at  n  =  0  produce  extensional  wave  responses  that  have  been  slowed  to  190  m/ 
sec.  The  reduction  in  extensional  wave  speed  occurs  because  the  ALCRYN  coating  contributes 
mass  with  a  negligible  contribution  to  the  stiffness  of  the  composite.  A  comparison  of  the 
optical  phase  sensitivity  to  radial  pressure  and  to  longitudinal  shear  stress  reveals  a  shear  stress 
response  50  to  70  dB  greater  than  the  radial  pressure  response.  This  result  is  consistent  with  that 
for  the  bare  fiber. 


The  Poisson  ratio  exerts  an  interesting  effect  on  the  dynamic  behavior  of  the  cylinder.  As  the 
value  of  the  Poisson  ratio  increases,  the  relative  difference  between  the  value  of  10Log(e^/P o) 
yt  =  0  and  the  value  of  XQhogitJPof  above  the  extensional  wave  resonance  (flat  plateau) 
increases  (figure  41).  In  the  limit  of  =  0,  the  extensional  wave  vanishes  under  a  P^ 
excitation.  This  behavior  is  consistent  with  the  fact  that  the  Pg  excitation  develops  an 
extensional  wave  response,  which  is  a  transverse  response  to  the  applied  force,  only  when  there  is 
a  means  of  coupling  energy  from  one  direction  into  the  corresponding  perpendicular  direction. 

Similar  to  the  effect  that  an  increasing  value  of  the  Poisson  ratio  has  on  the  single  cylinder,  a 
second  cylinder  with  differing  material  properties  from  the  first  serves  to  increase  the  strain  and 
optical  phase  sensitivity  in  comparison  to  the  bare  optical  fiber  in  the  region  of  the  extensional 
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wave  resonance  peak  (figure  83).  Very  little  change  occurs  near  A:  =  0  for  this  infinite  cylinder 
formulation. 
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APPENDIX  A 

BESSEL  FUNCTIONS  OF  COMPLEX  ARGUMENT,  SERIES 

REPRESENTATION 

The  inclusion  of  structural  damping  in  the  cylinder  necessitates  the  evaluation  of  Bessel 
functions  of  complex  argument.  The  series  representation  of  these  functions  is  then  required. 
From  Abramowitz  and  Stegun,^  equations  (A-1)  through  (A-6)  are  used  for  complex  argument  z, 
where  z  ranges  from  0  ^  z  ^  3 .  The  circumferential  order  number  n  appears  in  the  equations  as 
weU. 

The  Gamma^  function  is  given  by 

r(n  +  K+l)  =  in  +  K)\.  (A-1) 


The  Psi^  or  Digamma  function  is 

H  -  1 

v(i)  =  -y>  v(«)  =  -Y+  ^  . 

k=  1 

The  value  used  for  Euler’s  constant  is  y=  0.577215664901532860606512. 


The  Bessel  function  of  the  first^  kind  is 


(A-2) 


k  =  0 


((4)= 


/:!r(Az  +  ^+  1)  ' 


(A-3) 


1 .  M.  Abramowitz  and  L  A.  Stegun,  Handbook  of  Mathematical  Functions,  U.S.  Government  Printing  Office, 
Washington,  DC,  June  1964. 

2.  M.  Abramowitz  and  I.  A.  Stegun,  p.  255,  equation  6.1.6. 

3.  M.  Abramowitz  and  I.  A.  Stegun,  p.  258,  equation  6.3.2. 

4.  M.  Abramowitz  and  I.  A.  Stegun,  p.  360,  equation  9.1.10. 
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The  Bessel  function  of  the  second^  kind  over  the  range  0  <  z  ^  3  is 


J^k=o  ^ 


f-1"  “ 

k  =  0 


(A-4) 


The  Hankel  function^  of  the  first  kind  is  given  by 


(A-5) 


The  expression^  used  for  the  modified  Bessel  function  /„  over  the  range  0  ^  z  ^  3  is 


M  2r 

T  Z"^ 


[2)  E  )t!r(n  +  it+ 1)  ■ 


(A-6) 


The  expression^  used  for  the  modified  Bessel  function  K„  over  the  range  0  <  z  ^  3  is 


+  (-l)"l^2j  J  L  (V(^+l) +¥(«  +  ^+l)) 


lz2V 


(A-7) 


5.  M.  Abramowitz  and  I.  A.  Stegim,  p.  360,  equation  9,1.1 1. 

6.  M.  Abramowitz  and  I.  A.  Stegun,  p.  358,  equation  9.1.3. 

7.  M.  Abramowitz  and  I.  A.  Stegun,  p.  375,  equation  9.6.10. 

8.  M.  Abramowitz  and  I.  A.  Stegun,  p.  375,  equation  9.6.1 1. 
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9.  G.  A.  Kom  and  T.  M.  Korn,  Mathematical  Handbook  for  Scientists  and  Engineers,  McGraw-Hill  Book 
Company,  New  York,  pp.  868-869,  equations  21.8-44  and  21.8.45. 
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APPENDIX  B 


DERIVATIVES  OF  BESSEL  FUNCTIONS 


An  expression  for  the  first  derivative  of  the  Bessel  functions  needed  for  the  evaluation  of  the 
displacement  potentials  is  taken  from  Arpaci,^  and  is  listed  here  as  equations  (B-1)  and  (B-2): 


^(ZJax))  =-aZ^^^(ax)+\(ax),  Z  = 


(B-1) 


^(IJax))  =  a/„^i(ax)  +^/„(ax)  . 


(B-2) 


The  argument  of  the  Bessel  functions  is  composed  of  a  (a  constant)  and  x  (a  complex  number). 

An  expression  for  the  second  derivative  of  the  Bessel  functions  listed  in  equation  (B-1)  is 
derived  by  differentiating  this  equation  again  with  respect  to  x,  which  requires  the  use  of  the  prod¬ 
uct  rule  on  the  second  term  of  the  right-hand  side  of  equation  (B-1)  and  the  following  effect  of  the 
recursion  relationships  for  Bessel  functions: 


z  =  j,y,k,h^^\h^^^ 


.  (B-3) 


The  final  expression  for  the  second  derivative  is 


^ZJax)  =  ^(Z^_^{ax)  -  (2Z„(a;c)  +Z^^^iax))) 


z  =  j,y,k,h^^\h^^^ 


(B-4) 


1.  V.  S.  Arpaci,  Conduction  Heat  Transfer,  Addison-Wesley  Publishing  Company,  Reading,  MA,  p.  139, 
equation  (3-139). 
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